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I We study the generation and evolution of density perturbations and peculiar velocities due 
0\ ■ 

' to primordial magnetic fields. We assume that a random magnetic field was present before 
recombination and follow the field's effect on the baryon fluid starting at recombination. 
2 ■ We find that magnetic fields generate growing density perturbations on length scales larger 
^ I than the magnetic Jeans length, A^, and damped oscillations for scales smaller than A^. 
^ ■ For small wavenumbers k (large length scales), we find the magnetic field-induced density 
^ I power spectrum generally scales as /c^. We derive the magnetic Jeans length explicitly by 
■ - - ■ including the back-reaction of the velocity field onto the magnetic field by decomposing 
the magnetic field into a force- free background field and perturbations about it. Depending 
on the strength of the magnetic field and the ultraviolet cutoff of its spectrum, structure 
can be generated on small or intermediate scales early in the history of the universe. For 
a present rms magnetic field of 10"^'^ G on intergalactic scales, we find that perturbations 
on galactic scales could have gone non-linear at 2 ~ 6. Finally, we discuss how primordial 
magnetic fields affect scenarios of structure formation with non-baryonic dark matter. 
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1. Introduction 



The past decade has seen tremendous growth in our observational picture of the uni- 
verse. The Cosmic Background Explorer (Smoot et al. 1992) and other cosmic background 
experiments have shown that the large-scale clustering seen in galaxy surveys is consistent 
with a primordial origin for density perturbations. On the largest scales, where density 
perturbations are linear (i.e., rms variations in the density are smaller than the mean), mi- 
crowave anisotropy observations point toward an approximately Harrison-Zel'dovich spec- 
trum of initial density perturbations {P{k) ~ oc fc", with n ~ 1; Smoot et al. 1992; 
Ganga et al. 1993). Such a spectrum arises naturally in inflationary theories as well as in 
models based on topological defects. 

On scales smaller than ~ 8h~^ Mpc (where h is the Hubble constant in units of 100 
km/sec/Mpc), galaxy clustering is non-linear, and one needs to rely on numerical studies 
to get insight into the physics of cluster and galaxy formation. (Hereafter, we set h = 1 
unless otherwise noted.) Not only do the density perturbations become non-linear, but the 
complexity of the physics involved escalates as hydrodynamical effects become important. 
In this paper, we show that an element of this increased complexity that is often neglected, 
namely magnetic fields, may play a key role in the formation of structure in the non-linear 
regime. 

Interest in the possibility of primordial fields has also been rekindled by recent studies 
of galactic dynamos, which suggest that the "classical" dynamo mechanism for galaxies, 
i.e., "cK — a;" mean field dynamos, provide an inappropriate description of the origins of 
galactic magnetic fields (Ko & Parker 1989; Rosner & DeLuca 1989; Kulsrud & Anderson 
1992; Chakrabarti, Rosner, & Vainshtein 1994). A variety of difficulties have been clearly 
recognized: To begin with, virtually all extant galactic dynamo models are kinematic, and 
hence cannot properly describe the evolution of the /? (= Sirp/B'^, p the gas pressure) fa 1 
galactic magnetic field. Buoyancy is typically parametrized, and not explicitly taken into 
account (despite the fact that it is probably the dominant process by which magnetic fields 
generated in the solar interior and the galactic disk interior reach the outer "surfaces" of 
these objects; Parker 1979.). The fields are assumed to be smooth; that is, there is usually 
no account taken of the observed filamentation and other fine structure evident in both 
solar magnetic fields and magnetic fields in our galaxy (although in the galactic case, it 
is not certain that the structure is in the field rather than in, for example, the particle 
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density distribution). The galactic fluid is viewed as a single-component gas; this means 
that the effects of neutrals (including such phenomena as ambipolar diffusion) are ignored, 
contrary to what one would expect (in contrast, cf. Kulsrud 1988; Zweibel 1988). Finally, 
there is some considerable doubt as to the actual value of the turbulent (eddy) diffusivity: 
It has been argued by Ko & Parker (1989) that conditions appropriate to magnetic field 
generation and turbulent diffusion occur rather sporadically, and so are not the norm (also 
Parker 1992); and various authors have recently argued that the concept of turbulent 
diffusion for magnetic fields is itself fiawed (cf. Vainshtein & Rosner 1991; Kulsrud & 
Anderson 1992; Cattaneo & Vainshtein 1991; Cattaneo 1994; Gruzinov & Diamond 1994). 
Thus, although it is by no means established that dynamo models to explain galactic fields 
cannot be constructed (indeed, see Parker 1992 for a recent new proposal), it is by the 
same token no longer certain that extremely weak primordial "seed fields" are sufficient to 
account for the observed fields in galaxies. 

Of great relevance to understanding the role of magnetic fields in galaxy formation are 
observations of intergalactic magnetic fields and magnetic fields at high redshifts. Reports 
of Faraday rotation associated with high-redshift Lyman-a absorption systems (Kronberg 
& Perry 1982, Wolfe 1988; Wolfe, Lanzetta, & Oren 1991, Kronberg, Perry, & Zukowski, 
1992) suggest that dynamically significant magnetic fields may be present in condensations 
at high redshift. Together with observations of strong magnetic fields in clusters (Kronberg 
1994), these observations support the idea that magnetic fields play a dynamical role in 
the evolution of structure. 

The notion that magnetic fields may play an important role in structuring the universe 
is not new. The most detailed study was done by Wasserman (1978), who assumed the 
existence of a random magnetic field at recombination and, by treating it as a "source 
term," showed that it may act as a source for galaxy-scale density fluctuations. These 
calculations did not take into account the possibility of fluid back-reactions to the Lorentz 
force via the induction equation. Hence, the temporal evolution of the magnetic fields was 
entirely due to the Hubble expansion, and a magnetic Jeans length could not be derived. 
As we discuss below, Wasserman's approach is valid for length scales much larger than 
the magnetic Jeans length, Xb, and fiuid velocities much smaller than the Alfven speed, 
va = B/^/Airp. For velocities approaching va, the back-reaction of the fiuid onto the field 
cannot be neglected. We study the v va regime by assuming that the background field 
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is close to force free; this treatment aUows us to derive the magnetic Jeans length within 
a linear analysis. 

In order to describe the implications of random magnetic fields present at recombina- 
tion for structure formation, we consider the coupled evolution of density perturbations, 
peculiar velocities, and magnetic fields for the following two cases. First, we study the 
effect of a background magnetic field to first order in the velocity and density fields (fol- 
lowing Wasserman 1978) and deduce the spectrum of density perturbations generated on 
scales in which the back-reaction onto the field is small, i.e., on scales A » Ag. Second, 
we include the fluid back-reactions to the magnetic forces, and consider how this type 
of dynamics determines the power spectrum of the resulting velocity field and density 
fluctuations when the background magnetic field is close to force-free. Inclusion of the 
back-reaction onto the magnetic fleld allows us to derive the magnetic Jeans length for 
this problem, while the assumption of a force-free background field allows a linear treat- 
ment of scales A ~ A^. Thus, we present a consistent linear perturbation analysis of 
the combined magnetic-fiuid evolution equations (in the single-fluid approximation), and 
compute the present density fluctuations and vorticity under the assumption that random 
magnetic flelds existed before recombination. We show that the resulting spectrum for 
density perturbations has a general form which is insensitive to the magnetic fleld spectral 
index; on large scales the spectrum of density perturbations is too steep {P{k) oc k'^) to 
flt the observed spectrum, while on small scales magnetic flelds introduce a peak in the 
spectrum around /speak — iiiin(A;B, A;inax)) where = 2tt/\b and fcmax is the ultraviolet 
cutoff of the magnetic fleld spectrum. 

The outline of our paper is as follows: We present the basic magnetohydro dynamic 
equations used in our analysis, discuss plausible initial conditions, and carry out our lin- 
earization, in §2. In §3, we discuss the case studied by Wasserman and calculate the 
spectrum of density perturbations for this case. In §4, we consider the case of close to 
force-free background magnetic flelds and compute the time evolution and spectrum for 
the compressible modes, the resulting density fluctuations, and the incompressible modes. 
Our results are discussed in §5. For the sake of clarity, we have placed details of our 
analysis in the Appendices. 

2. The Perturbation Analysis 
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In this section, we review the physical conditions at the time of recombination, discuss 
the basic equations used in our analysis, and develop our perturbation scheme. 

We assume that random magnetic fields were present at the epoch of recombination 
and that these were formed through pre-recombination processes (e.g., Hogan 1983; Turner 
& Widrow 1988; Quashnock, Loeb, & Spergel 1989; Vaschaspati 1991; Ratra 1992a,b; 
Cheng & Olinto 1994). As the universe cooled through recombination, baryons decoupled 
from the background radiation, and the baryon Jeans length decreased from scales compa- 
rable to the Hubble scale (~ 100 Mpc) to ~ 10 kpc in comoving units. (Throughout this 
paper, we use A to describe comoving length scales and set the scale factor today to unity, 
R{to) = 1; physical length scales, I, at any time can be found by multiplying the comoving 
scale by the scale factor, l{t) = XR{t).) After recombination, baryons are free to condense 
on scales larger than the Jeans length, and will do so via gravitational instabilities if there 
are initial perturbations in the density field. Concurrently, magnetic fields that were frozen 
into the baryon-photon plasma before recombination will tend to relax into less tangled 
configurations once the baryons they are coupled to decouple from the photon background. 
Consequently, density perturbations in the baryons can be generated through the Lorentz 
force even if the density field is initially smooth, and the initial peculiar velocity field van- 
ishes at recombination. To understand the efi'ect of magnetic fields on the origin of density 
perturbations, we assume that no initial density perturbations or peculiar velocities were 
present at recombination, so that all subsequent density perturbations or velocities are 
induced by magnetic fields alone. (We address the more general case of combining initial 
density perturbations and magnetic field effects in a subsequent paper.) 

To follow the evolution of the density, peculiar velocity, and magnetic field after recom- 
bination, we write the basic one-fiuid magnetohydrodynamic (MHD) equations in comoving 
coordinates. 




(1) 



dtp + 3—p + 
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a.(«.B) = Ziiil^ (5) 

-ft 

(cf. Wasserman 1978), where ip is the gravitational potential, R is the scale factor, pf, = 
Pb{t) is the uniform background density, and all other symbols have their usual meaning 
(we set c = 1 throughout the manuscript). We have neglected all viscous and diffusive 
terms because the relevant Reynolds numbers are very large. 

2.1 Initial Conditions and Background Evolution 

We begin with the basic assumption that all baryonic matter is uniformly distributed 
at recombination (tree); with density Pb{trec) cind that this matter has zero peculiar velocity, 
v(x, tree) = 0. Furthermore, we assume that there is a magnetic field already present, 
B(x, tree) = Bi-ec(x), presumably created well before recombination; we posit that this 
magnetic field is randomly oriented on spatial scales smaller than the Hubble radius at 
recombination and has no ensemble averaged mean components on the Hubble scale; thus, 
we assume that (Brec(x)) = 0, while (|Brec(x)p) ^ 0, where the angular brackets mean 
the ensemble average. 

The above magnetic field does not significantly perturb the baryonic matter until 
photons and baryons decouple; once decoupling occurs, the unbalanced Lorentz forces act 
to disturb the smooth background density, Pb{t), leading to both density perturbations, 
5p{x,t), and peculiar velocities, v(x, t), of the baryonic fluid (Wasserman 1978). We 
decompose the fluid density, p(x, t), and the magnetic field, B(x, t) by writing 

p(x,t) = pfc(t) + 5p(x,t) = P6(t)[l + 5(x,t)], 
B{x,t) = Bb{x,t) + 6B{x,t) . 

Here B is the total field, B^ is the background random magnetic field, with 

R^ {tree) 



Bb{x,t) = Brec(x)- 



i?2(t) ' 



and 6B is the diff'crcnce between the total field and the background; in other words, Bi,{t) is 
simply the initial random field evolved only by the Hubble fiow, while SB is the additional 
field which results as the Lorentz force perturbs the baryonic fluid, and the fluid reacts 
back. 
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With these definitions, we assume that at t. 



S{x,trec) =0, 
v(x, tree) = , 

5B(x,W) =0. 

In addition, the temporal behavior of the background quantities is determined by the 
cosmological model we assume. To isolate the efi^ect of magnetic fields from other sources 
of density perturbation and to keep the analysis simple, we chose to study first the case 
of a fiat universe with a critical density of baryons (fibaryon = !)• In §3, §4, and §5, we 
discuss the effect on our results of including non-baryonic dark matter as the dominant 
component of the universe. 

For an Einstein-De Sitter model, the following relations hold during the time of in- 
terest (0 < ^ < Zrec ^ 1100): 

Pb{t)RHt) = Pb{to) = po, (6a) 
6TTGpb{t)t'^ = 1, (66) 
R{t) = (^/^o)'/^ (6c) 

where to is the age of the universe, and quantities with subscript are evaluated at the 
present time (we chose R{to) = 1). The background magnetic field satisfies 

Bb(x, t)R\t) = Brec(x)i?2(trec) = Bfa(x, to) = Bo(x) . {6d) 

It is also useful to define 

{Bl)^{\Bo{^)\') = 8nR{tfpB{t) , (7) 

where psit) is the average magnetic field energy density present in the universe at recom- 
bination redshifted to time t. 

In analogy with the ordinary Jeans length, Aj = Vs ^J'n j pG^ one can expect the 
magnetic Jeans length to be 

\b ~ va^Jt^IpG 
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(Peebles 1980). In §4, we show that the magnetic Jeans length, defined by the wavenumber 
where the transition between growing (or decaying) modes and oscillatory modes occurs, 
can be written as: 

The magnetic Jeans length written above can be constrained by requiring that, at the 
time of nucleosynthesis, the average energy density in the magnetic field, 

(|B(x,t)p) _ {Bl) 
P^^^^- 8^^ -8^^' 

be significantly less than the energy density in radiation. 



Pr{t) 
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where (7* is the number of relativistic degrees of freedom, and T is the temperature of the 
background radiation (Barrow 1976, see also Cheng, Schramm, & Truran 1994, Grasso & 
Rubinstein 1995, and Kernan, Starkman, & Vachaspati 1995). Since both energy densities 
redshift as it suffices to require that ps ^ Pr today, which gives: 

{B^)^ <4x 10"^ G. 

By assuming a critical density in baryons, po — 2x 10~^^ g cm""^, in eq. (8), we can write: 

Therefore, the constraint ps ^ Pr implies that the magnetic Jeans length today must 
satisfy lB{to) = Xb ^ 100 Mpc, i.e., 

(Note that if we had set Obaryon < 1 in the density used in eq. (8), the effect would be to 
increase the maximum A^, thus weakening the constraint.) 

It is interesting to note that, unlike the ordinary Jeans length, does not change 
with time; for t^ec < ^ < ^O; Isit) = XsRit)- 

There are several observational constraints on a primordial magnetic field, the most 
stringent being an upper limit on the Faraday rotation of light from distant quasars which 



requires that the average magnetic field on scales comparable to the present horizon, 
Bo{Hq^) ^ 10-11 G (Kronberg & Simard-Normandin 1976), where Hq^ ~ 3 Gpc. Using 
similar techniques to test the magnetic field on Mpc scales gives limits of So(Mpc) ^ 10"^ 
G, which is consistent with the limit one gets by assuming that the observed galactic 
magnetic field is solely due to a primordial field enhanced by the collapse of the Galaxy. 
Expanding the volume occupied by a galactic field of Bgai — 3//G in order to decrease the 
density from that in the Galaxy to the critical density gives an average field in the universe 
on the scale Ig ~ Mpc (the comoving scale that collapsed to form a galaxy) of Bo{lg) < 10~^ 
G (using 2 X IQ-^^g cm"^ for the average gas density in the Galaxy, 2 x 10~^^g cm~^ for 
the average gas density in the universe, and the assumption that the field is frozen in as 
the gas contracts, B oc p^/^). This argument is based on the simple assumption that the 
magnetic field is frozen-in and grows with the collapse as S oc p^/^. Kulsrud (1995) has 
argued that by taking into account the formation and rotation of the Galactic disk a seed 
field above -Bo(^g) ~ 10~^^ G may give rise to stronger galactic fields than observed. The 
dynamics of the Galactic field once fields close to or even above equipartition are reached 
has not been fully understood yet, so we will take the limit on Mpc scale fields to be 
So (Mpc) ^ 10"^ G below. 

These limits cannot be unambiguously translated into a limit on (Bq)^ since -Bo(/) 
refers to the average field on a particular scale I, for example, averaged with a window 
function (cf. eq. (42) below). This averaging procedure depends not only on the integrated 
power spectrum (Bq)^, but on the functional form of the power spectrum B^{k) as well, 
unless {Bq)^ ^ 10~^ G. We return to these constraints during our discussion of magnetic 
field spectra in §3 and §4. 

2.2 The Perturbation Scheme 

The next step is to identify our "small" quantities, which will fix the ordering of the 
perturbation scheme. In the spirit of a linearized theory, we shall assume that the density 
perturbations resulting from the Lorentz force are small, i.e., that 

(5< 1. 

Similarly, we assume that the induced peculiar velocities and magnetic fields are small, 
e.g., we assume that 

^^rmsT ^ , , \SB\ 

— - — -C 1 , and -C 1 , 



where r is the characteristic time scale of the flow, v^ms is the rms flow speed, and L is the 
characteristic length scale of the flow at time t^ec + T. Using these scaling relationships, we 
linearize eqs. (l)-(5), noting that since we are primarily interested in wavelengths larger 
than the ordinary Jeans length, the pressure term may be ignored (Peebles 1980; 1993). 
Then, upon retaining leading order terms, we obtain 

Vy; , (V X Bb) X Bb 

5t5+^ = 0, (12) 

= ^TTR^GpbS, (13) 

V • Bfe = V • 5B = 0, (14) 

dtiR^Bb) = 0, (15) 

Jrt 

Note that eqs. (15) and (16) allow us to clearly distinguish between the background field 
Bfo and the perturbed field 5B. 

The evolution of v due to the Lorentz force depends non- linearly on B. If we assume 
that the background field at recombination is not force-free, (V x B5) x B5 7^ 0, the 
background field will generate a velocity field with ^^v oc (V x B5) x B5. The fluid motion 
will affect the background field through cq. (16), but to higher order in perturbation 
theory. The fluid back reaction can be neglected to lowest order as long as |5B| ^ IB5I. 
This is the case studied by Wasserman (1978), who derived the time evolution of density 
perturbations, which we briefly review in §3.1. In §3.2 we derive the spectrum of density 
perturbations for this case. 

The effect of the fluid motions on the magnetic fleld cannot be neglected when the 
fluid velocities become comparable to the Alfven speed, va = B/y/47rp. This can be seen 
by studying the effect of one A;-mode of B^ in eqs. (11) and (16). Eq. (11) gives 

V kt 
Va R 



while eq. (16) implies 



SB kt V f V 

Bh R VA \VA 
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Therefore, 55 <^ if v <^ v^, and the back-reaction of the fluid onto the field can be 
neglected for scales A » A^. 

This first order description in which the background field is frozen-in (only redshifts) 
is physically reasonable for scales As -C A ^ A^, in which a tangled background field 
"recently" entered the horizon and had no chance to relax to a force-free configuration. 
For smaller scales, A ~ As, an alternative description is necessary. At recombination 
Xh is much larger than A^, for magnetic fields of interest ((-Bq)^ < 10 ~^ G implies 
Ab < 25 kpc), so scales close to A^ have been inside the horizon for some time and 
will tend to relax to force-free configurations.^ At recombination, intermediate scales 
(As <^ A <^ Ah) will be relaxing toward a force- free configuration while scales close 
to Ab behave like oscillatory modes about a force-free configuration. For this range of 
scales, a more appropriate description of the evolution of the fiuid and magnetic field 
can be constructed by assuming that the background field has a force-free component 
on intermediate to large scales, B/j-, and small scale perturbations about this force-free 
configuration, b, 

B, = Bff + h. (17) 

By assuming that |B//| >> |b|, we can linearize eqs. (1) and (5) by keeping terms to first 
order in v and b, and write 

Vy^ , (VxB^^)xb+(Vxb)xB^y 

dt{R^Bff) = 0, (19) 
^^^^,^^^V><(vx^_ (20) 

Jri 



^ There is a substantial literature on the subject of how a magnetized plasma might 
relax in the presence of constraints, starting with the classic papers of Woltjer (1958) and 
Taylor (1974). For low-/? plasmas, Taylor argued that relaxation occurs such that the total 
helicity, H = J dr^ A ■ B (with A the vector potential corresponding to the magnetic field 
B) is constant; in that case, such a system will relax to its lowest energy state consistent 
with this fixed value of H, with the important property that the magnetic field is force- 
free, i.e., V X S = X{'ip)B, where ip is the fiiix function corresponding to B. Thus, as long 
as one can ignore the dynamical effects of the baryons, one would expect the fields to relax 
to a force-free state. 



11 



Eqs. (18)-(20) together with eqs. (12)-(14) allow us to study the behavior of the 
coupled equations when the fluid back-reaction becomes significant and the magnetic Jeans 
length can be derived. This approach is valid when a separation of scales between the 
force-free component and the perturbations about it is possible. To get the spectrum of 
the generated density perturbations in this case, we need to specify both the spectrum of 
Bff and b at recombination (B(,(x, = B/j-(x, tree) + b(x, t^ec) and (5B(x, t) = 

b(x, t) — b(x, tree))- Following this scheme, we deduce the time evolution and find Xb in 
§4.1, and discuss the spectrum of the generated density perturbations in §4.2. 

The physical picture and plan of this paper is as follows: Starting at recombination, 
we follow the effect of a random magnetic field present at recombination on the velocity 
and density fields. For scales Ab <S A ^ A//, we use the linear approach developed by 
Wasserman, who derived the time evolution of these fields. We recapitulate Wasserman's 
results in §3.1 and derive in §3.2 the spectrum of generated density perturbations by 
assuming either a power law or a delta function spectrum for the magnetic field. We show 
that the power spectrum of density perturbations is too steep to fit the observed clustering 
of galaxies on large scales. For A ~ Ab, we derive in §4 the time evolution and spectra for 
density and velocity fields by assuming a force-free background field and either a power 
law or delta function spectrum for the magnetic field perturbations about the force-free 
background. 

3. Fixed Background Magnetic Field 

In this section, we first discuss the time evolution of the compressible modes for the 
case in which the field has not relaxed to a force-free background configuration and acts 
as a fixed comoving source of density perturbations that only redshifts in time. This 
approach is valid to first order in perturbation theory and for scales X ^ Xb- We then 
solve for the spectrum of generated velocity field and density perturbations and discuss 
the implications. 

3.1 Time Evolution 

When the back-reaction of the fluid onto the background magnetic field may be ne- 
glected, (i.e., for A » A^), we can solve for the time evolution of 5{x,t) by taking the 
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divergence of eq. (11) and using eqs. (12) and (13) (Wasserman 1978). This gives: 



^ ^^o^^A A n X V • [(V X Bb) X Bj,] 
du5 + 2-d,5 - A.Gp,5 = ^—^^ . 

Making use of the time dependence of R, pb, and B;,, through eq. (6), we can write: 

2 ^toV V-[(VxBo)xBo] 
^«'^+3^^*^-3^^=UJ ^0 • 



(21) 



(22) 



(Note that in our convention po = Pb{to) while Wasserman uses po{t) for the background 
density.) 

The time evolution of 5 is then given by 



5(x,t) = t^ 



2V • [(V X Bo) X Bo] 



47rpo 



9 / t 

10 [tZ, 



'^ 3 tree 3 

5~r ~ 2 



(23) 



The time evolution in this case is separable from the spatial dependence such that every 
A;-mode varies the same way in time; growing modes scale as t^/^ or as R{t), while decaying 
modes scale as t~^. 

The time evolution of the compressible modes can be found via eq. (8), namely 
3WV • [(V X Bo) X Bo] / to 



V-v(x,t) = 



free 



207rpo V^rec 

Finally, the time evolution of the incompressible mode is given by 

3trecV X [(V X Bo) X Bo] / to 



(24) 



V X v(x, t) = 



A-KPQ 



tr 



f \ 3 

T 



2 

f ^ 3 



(25) 



Unlike compressible modes, incompressible modes have no growing modes during the mat- 
ter dominated era. 

3.2 Spectra of Density Perturbations and Velocity Fields 

In order to obtain the spectral dependence of 6 and v for a given spectrum of magnetic 
fields at recombination, we define the Fourier-transforms of the fiuid variables in comoving 
coordinates, 

5(x, t) = J exp(zk • x)5(k, t) , (26) 
v(x, t) = J exp(zk • x)v(k, t) , (27) 
Bo(x) = /d^k exp(ik-x)B(k). (28) 
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Brec(x) is assumed to be homogeneous and isotropic (and, obviously, so is Bo(x)); therefore 
B(k) obeys the relation (cf. Kraichnan & Nagarajan 1967; Moffatt 1978) 



(5,(ki)5;(k2)) = S\k, - k,) { S,, - ] , (29) 



kl ) 2 



where i and j label the z-th and j-th components of the vector B(k). Note that 

{BD = 47r y dk^B^k) = G {^-^^ ' • (30) 

For the case of a fixed magnetic background field considered in this section, the spectral 
dependence of 5 is solely due to the initial conditions at recombination. In §4, we discuss the 
case in which back-reactions become important and the time dependence is a function of the 
wavenumber. Here, it is suflBcient to determine the spectral shape of dtt5 at recombination 
which is given by eqs. (12) and (24), 

V-9,v(x,t,ec) V ■ [(V X Bo) X Bo] 
-Krec 47rpoi?^ec 

We would like to determine the spectral dependence of the ensemble average of 5 
by assuming a spectrum B'^{k) for a random magnetic field at recombination as in eq. 
(29). The ensemble average of both sides of eq. (31) vanish, since there are as many 
positive and negative variations of du5., but the relevant quantities to be computed are the 
ensemble averages of the square of the fluid quantities. The root mean square (rms) of the 
Fourier-transform of the density field is defined by 



^rms(/c,t)^7(|5(k,t)|2). 



Since, in the fixed background field case, the time and spectral dependencies are separable, 
we can write 

Srtns{k,t) = c{k)T{t) , 

where 

and c{k) is determined by eq. (31). The spectral dependence of 5rms is related to the rms 
of the curl-free component of the velocity field which we define by 



^^||rms(^,^) = ^(l^||(k, t)|2) , 
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where k- v(k, t) = kv\\ (k, t), and the subscript || represents the compressible (longitudinal) 
component of v. Similarly, we define 



V±rins{k,t) = , 

where k x v(k, t) = kv±{k,t), and the subscript _L represents the incompressible (trans- 
verse) component of v. 

The detailed derivation of c{k) and 'i'||rms(^) ^) as a function of B^{k) is provided in 
Appendix A; here we only write the solution, 

= ^^||rms(fe, W) k , (33) 



and 



Vnrms{k,trecf k^ = { ] JT^ dki dfi ^' ^'^ 



{2A;^A;?// + A;^A;t(l - S/x^) + 2A;^/c^^^} , (34) 

where F is a volume factor and k • ki = kki/j,. 

Since a given /c-mode of the velocity field is a non-linear convolution of the magnetic 
field spectrum, different A;-modes of B{k) will contribute to a given mode of v{k). In 
particular, if the spectrum for the magnetic field has an ultraviolet cutoff, /cmax, the velocity 
field modes will be non-zero for k < 2kma,x- In principle, the lowest /c-mode excited through 
eq.(31) is A; = 0, even if the magnetic field spectrum has an infrared cutoff, but, in practice, 
the Hubble radius at each time will provide an effective infrared cutoff for the density 
perturbations excited by the magnetic field. As we show below, the power in low A;-modes 
is too small for the infrared cutoff to be relevant. 

In what follows, we discuss two Ansatze for the functional form of the magnetic field 
spectrum: power laws and delta functions. The delta function Ansatz is simple to calculate 
and can be used to relate our results to those of Wasserman (1978) and to the more realistic 
case of the power law Ansatz. 

3.2.1 Delta Function Spectra for B'^{k) 

If the magnetic field spectrum at recombination is a delta function, we can write 

BHk)^{Bl)'-^i^. (35) 
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The initial acceleration becomes for k < 2k^, 

= ^(sSt^ (36) 

and zero otherwise (see Appendix A). Using eq. (33) and ks from eq. (30), we find for the 
density spectrum 

25VV ki 

^^^^= 384^^- ^'^^ 
We can now calculate the power spectrum of density perturbations, P{k,t), and the 
variance {rms power per logarithmic wavenumber interval), A{k,t), using the following 
definitions 

(27r)3 



P{k,t)^'-^{\6{k,tr) (38) 



such that. 



and 



such that. 



(|5(x,t)|2)^ I d'kPik,t), 



327r^ 

A{k,tf = 4nk^P{k,t) = ^/c^(|(J(/c,t)|2) (39) 



(|5(x,t)|2) = J A{k,tfd\nk. 
For the delta function magnetic field spectrum, we obtain 

^ „ , 25V2 f kV ks , ^ 

The resulting power spectrum of magnetic field-generated density perturbations can 
be compared to observations of galaxy clustering, if we assume that galaxies trace the mass 
density distribution. The power spectrum of galaxy clustering has been measured over the 
range 0.1 Mpc < 2TTk-^ < 10^ Mpc (e.g., Geller & Huchra 1989; Efstathiou et al. 1990; 
Maddox et al. 1990; Collins, Nichol, & Lumsden 1992; Fisher et al. 1993), while information 
from anisotropies in the background radiation reaches scales comparable to the present 
Hubble radius 27TkQ^ ~ 3 x 10^ Mpc (Smoot et al. 1992; Ganga, Cheng, Meyer, & Page 
1993). In the absence of bias, an estimate of the present-day scale of non-linear clustering 
can be made by estimating the scale at which the rms galaxy fiuctuations are unity; for 
optically selected galaxies this is lni{to) = A„/ ~ 8 Mpc. On large scales (O.OlMpc"^ < 
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k < kni), the observed galaxy power spectrum is consistent with a power law, P{k) oc A;", 
n ~ — 1 with a hint of a bend to larger n at the larger scales. COBE suggests that 
P{k) (X k on the largest scales. Therefore, magnetic field-induced perturbations (with a 
delta function magnetic field spectrum) have too steep a spectrum (Psik) oc k^) to agree 
with observations on large scales. As we discuss below, a similar behavior is found if 
B'^(k) is a power law (in that case, P{k) oc k'^)] therefore, magnetic field- induced density 
perturbations cannot reproduce the observations of structure on very large scales. 

On smaller scales (k > 27r /Xni), magnetic field-induced perturbations may play an 
important role. In particular, they are of interest if B^{k) is such that As{k,to) > 1 
for cosmologically relevant scales, say between clusters of galaxies {kd ^ 27r/2 Mpc) and 
globular clusters {kgi ~ 27r/10~^ Mpc). As As{k,t) approaches 1, the linear treatment 
used above breaks down. However, we can make use of our linear solution to approxi- 
mately estimate the epoch, tni{k), that a particular scale k becomes non- linear by setting 



The first scale to go non-linear in this case is the smallest wavelength allowed or largest 
wavenumber k — 2k^. The time at which perturbations with k = 2k^ become non-linear 
is then tni{2k^) ~ 1.8 tree {kB/2k^Y . Requiring that tni{2k^) < to implies 2k^ ^ 3.7 x 
10~^/cb. For example, suppose = O.lkB, which gives tni{2k^) ~ 1.8 x lO^trec or 
redshift Zni{2k^) ~ 6. In this case, choosing the galaxy scale to go non-linear at 2; ~ 6, 
we get 2A;* ~ 27r/ Mpc, ks — 27r/0.1 Mpc, and the constraint eq. (10) is satisfied with 
(Bq)-*^/^ ~ 4 X 10~^ G. This scenario would correspond to the formation of galaxies around 
redshift 6. If instead, we choose cluster scales to go non- linear at redshift Zni = 1, then 
2fc* ~ 27r/2 Mpc and ks — 27r/0.1 Mpc, the same as the previous example. A final example 
would be to choose the non-linear scale to be at present, tni{2k^) = to and 2k^ — 27r/8 
Mpc, then ks = 2n/0.3 Mpc, which corresponds to a field {B^y/'^ ^ 10~^ G. 

The examples given above satisfy the constraint that Pb P-y- As we mentioned 
before, there are other constraints on the contribution of B'^{k, t) to a given scale in which 
Faraday rotation measurements can be made. Unfortunately, these measurements can only 
be translated to a limit on magnetic fields when assumptions about the density of ionized 



A{k,tnl{k)) ^ 1. 

Setting As{k,tni{k)) ~ 1, we find that tniik) satisfies 




(41) 



17 



gas in the line of site are made. Since we have assumed a critical density of baryons, 
applying the constraints from these measurements can be somewhat ambiguous in the 
present discussion. When we address similar effects in the presence of dark matter some of 
the ambiguity will not be present (e.g., the dependence on Obaryon can be kept explicitly). 

Keeping in mind these limitations, we briefly discuss the implications of these con- 
straints for the examples given above. To apply these observational constraints, we define 
the average field on scale I by 



where W{'x. — x.') is a window function which smoothes the magnetic field on scale I and 
W{kl) is its Fourier transform. Using a Gaussian window function, VF(x — x') = exp(— |x — 



For the three examples discussed above, the formation of galaxies at 2; = 6 and the forma- 
tion of clusters at 2; = 1 satisfy eq. (43), while the parameters chosen for the formation 
of structure on scales lni{to) today do not satisfy this constraint. 

The examples discussed above demonstrate how non-linear structures within a limited 
but relevant range of scales can be formed at reasonable redshifts if magnetic fields satis- 
fying observational limits were present at recombination. Because of the steep spectrum, 
the scales influenced by magnetic flelds are primarily in the non-linear regime, ultimately 
requiring a detailed numerical study. Although the Ansatz used above for the magnetic 
field spectrum is not realistic, some of the results obtained above are very similar to the 
power-law spectrum discussed below. 

Before leaving this section, we note that Wasserman (1978) discussed the case in which 
B'^{k) is sharply peaked around k = 27r/xG (his xq corresponds to our Ig) and wrote (his 




(42) 



x'|V2Z2) in eq. (36), we find 




Requiring that B{1, to) < 10 ^ G for Z = ~ Mpc, we can write 




(43) 
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eq. (26)) 

- 2 

(|V • [(V X B,ec) X Brec]|')^ ^ ^ (^) 

as an estimate for the effect of magnetic fields. His result corresponds to an average over 
the volume 

(IV ■ [(V X B.ec) X Brecin = ^'^'^ ^'^^'^^ J d'^{\V ■ v(x, t)\') 

which, through Parseval's theorem, can be re-written as 

(IV ■ [(V X Brec) X Brec]|') = illlM^^l^ (27r)3 | d'k{\h ■ ^(k,t)p) . 

Choosing — 2t\ jxQ in the delta function Ansatz for B^{k) and using eq. (36), we recover 
Wasserman's eq. (26). 

3.2.2 Power Law Spectra for B'^{k) 

A variety of mechanisms for generating magnetic fields before recombination have 
been proposed (Hogan 1983; Turner & Widrow 1988; Quashnock, Loeb, & Spergel 1989; 
Vaschaspati 1991; Ratra 1992a,b; Dolgov & Silk 1993, Dolgov 1993; Cheng & Olinto 
1994), but consensus on a well-motivated scenario is still lacking. In general, most models 
generate a power law spectrum, B'^{k) oc /c^, with a cutoff on small scales (a typical 
cutoff scale is some fraction of the Hubble radius when the field was generated). Of 
particular interest is the case of a white noise {q = 0) spectrum, which would result from 
magnetic fields generated with similar strengths but random directions within each Hubble 
volume during a phase transition in the early universe (e.g., Hogan 1983). After the phase 
transition, this can be viewed as a random walk of field lines, with stepsize of the order of 
the coherence length (the Hubble radius at the phase transition, or some fraction of it). 

The magnetic field spectrum generated in the early universe will evolve differently on 
different scales. On very large scales, the field is frozen into the fluid and only redshifts 
with the expansion of the universe, B oc R-^. On very small scales, the finite plasma con- 
ductivity allows diffusion of the fleld within the plasma. On intermediate scales, damping 
of the magnetohydrodynamic modes as the neutrinos decouple and later as the photons 
decouple will change the effective cutoff of the magnetic field spectrum to scales as large 
as the Silk scale (Jedamzik, Katalinic, & Olinto 1995). 
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Here, we restrict our attention to the evolution starting at the end of recombination 
for a magnetic field spectrum parametrized by the power law index q and an ultraviolet 
cutoff fcmax- We assumc that, for kmin < ^ < 



"max) 



B\k) = Ak'i, (44) 

where ^ is a constant, ^ ~ (q+S) (Bq) / (for small fcmin). To calculate ■U||i.ms(^5 W) an- 
alytically, we make the further assumption that the magnetic field is Gaussian distributed. 

The details of the tedious algebra are left to Appendix A. To leading order in /c//Cmax <^ 1 
and for integer spectral index q between -1 and 6, we obtain the generic result that 

i\\rms{k, tree) - ^ \ — , (45) 



"'max 



where 

+ 3) 



e„ = 



2Vl5(2g + 3) ■ 
Using eqs. (33) and (45), we find 

25^^ eoA;2 



c{k) 



"'^'I'max 



I 3 



^q%t)c.'^eA^\ \-^\ Tit) , (46) 



48 " \kB 
and 

The spectrum of generated density perturbations is almost independent of the mag- 
netic field spectral index: ~ 1, for — 1 < < 6. (This range in q spans most of the 
proposed primordial magnetic field spectra, with the exception of Ratra (1992a,b) who 
proposed indices up to g = —3; wc return to this issue in §4.2). In contrast, there is a 
strong dependence on the ultraviolet cutoff /Cmax, which plays a role similar to /c* in the 
delta function Ansatz. The amplitude of the spectrum of perturbations is determined by 
kB and A;niax- For the power-law indices discussed above, the dependence on an infrared 
cutoff, Ajmin, is negligible unless A;inin is larger than the wavenumber of interest. 
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The power spectrum is a power law, Pq oc k'^, for scales inside the horizon at recombi- 
nation, i.e., for k > kj-ec = STrffreo is steeper for k < k^ec, with a cutoff at 2A;niax- The 
peak power is around k ~ 2A;inax although the exact behavior cannot be obtained from eq. 
(45), since we neglected terms of higher order in fc/fcmax- 

An interesting property of the spectrum Pg is that, unlike the case of the delta function 
Ansatz, if we evaluate it at /cb, Pqiks) does not depend on ks- This implies that if we 
extrapolate our solution to the limit k ^ ks, Pq would have a fixed shape and amplitude 
as /cs is changed, only shifting horizontally on a Pq vs. k plot. 

Before exploring the relevant ranges in A;niax and ks for the formation of structure, 
we discuss the observational constraints on both parameters. The constraint on ks in eq. 
(10) is unchanged, since it is independent of the magnetic field power spectrum, while the 
constraint on kma.x depends on the spectral index q. For q= 0, the average field on scales 
I is given by (from eqs. (42) and (44)), 



S(Mo) = ,/^^, (48) 



"max" 



where we use the same Gaussian window function as in the previous section to smooth the 
field on scales I. Again, for I ~ Ig, B{lg,to) ^ 10~^ G and we get 



which is more restrictive than eq. (43). As q increases, the constraint becomes less severe, 
since the steeper the spectrum the less it contributes to Bi. 

We now return to the ranges in fcmax and ks that are relevant to structure formation. 
In an analogous way to the delta function case, small scales reach non-linear variance 
(A > 1) earlier than large scales, and the first scales to become non-linear have k ~ 2/cinax- 
We again define tni{k), such that Aq{k,tni) — 1, which implies: 



2 3 7 

tni{k)Y 1.5 fk 

max \ / i^B 



tree J \ ks J V ^ 



(50) 



Requiring that tni{k ~ 2/cinax) ^ leads to the constraint /cmax ^ 10~^/cb- This is 
an approximate estimate for the lower limit for A;max, since eq. (45) was derived in the 
limit of A; <^ A;niax- If we again choose 2A;niax = 27r/ Mpc and require that this scale go 
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non-linear at redshift Zni — 6, we get ks = 27r/0.11 Mpc which satisfies both constraints 
from eq. (10) and eq. (49). This would correspond to the formation of galaxies at redshift 
6 for a background field (B^) ~ 4.4 x IQ-^ G and a B{lg) just below 10"^ G. If instead, 
we require that magnetic fields form clusters of galaxies at redshift Zni = 1, that would 
correspond to 2/cjnax = 27r/2 Mpc and it would require ks = 27r/0.12 Mpc. This scenario 
satisfies eq. (10) easily and may or may not satisfy eq. (49) depending on the precise 
value for the scale at which the constraint applies and the value of the Hubble constant. 
(For example, if we reintroduce the appropriate /t-dependences, eq. (49) is satisfied for 
B < IQ-^ G and /i = 0.5 with Ig > 1.2 Mpc , or if /i = 0.8 then Ig > 1.7 Mpc.) A scenario 
for generating structure on scales which are going non-linear today would correspond to 
choosing the scale 2A;niax = 27r/8 Mpc to go non-linear at Zni = which satisfies eq. (10) 
but does not satisfy eq. (49). 

We see that with a more realistic choice for the power spectrum of magnetic fields at 
recombination, non-linear structures in the cosmologically relevant range of scales can be 
formed at reasonable redshifts for field strengths which are below the present observational 
limits. Some of the examples above have field strengths that are within the range of future 
observational capabilities (Kronberg 1994; Lee, Olinto, & Sigl 1995). 

We concentrated on galaxy and cluster scales in our examples above, but another 
possible consequence of magnetic fields at recombination is the formation of smaller objects, 
such as QSO's or Pop III stars, at very early times. The challenge in the formation of 
smaller objects is the survival of power on those scales through recombination, since photon 
drag damps magnetohydrodynamic modes up to the Silk mass (Jedamzik, Katalinic, & 
Olinto 1995). 

The scales on which magnetic fields generate structure are primarily in the non-linear 
regime today, which limits our ability to make precise predictions within linear perturbation 
theory. In the case of the baryonic universe studied above, the range of scales affected by 
magnetic fields is quite narrow, which suggests that objects formed via magnetic field 
perturbations may be biased with respect to large scale structures formed by primordial 
perturbations. 

For completeness, we derive the spectrum of generated incompressible modes. The 



22 



initial value of 5_Lrms(^) ^rec) can be obtained from eq. (6) by taking the curl, 

V X [(V X Bo) X Bo] 



VXV(X, tree) 



47rpo-Rre 



Again, as in the compressible case {v ±{'k, tree)) = 0. 

If we assume a power law Ansatz for the magnetic field spectrum, as in eq. (44), we 
find for the incompressible velocity field spectrum: 



. {q + 3){Bi)k I 28V 



(47r)3poi2L V 15(2g + 3)A;3^,/ 

Comparing to eq. (45), we see that the incompressible mode is initially excited to almost 
the same extent as the compressible mode. Although the two modes have comparable 
"initial" amplitude and spectrum, they have quite different temporal evolution and remain 
decoupled in the linear regime. 

4. Force-Free Background Magnetic Fields 

In this section, we consider an alternative perturbation scheme more appropriate for 
scales close to As. We study the effect of magnetic fields described by small scale pertur- 
bations about a large scale force-free background field. We derive the time evolution of 

the compressible and incompressible modes in §4.1 and in §4.2 we discuss the spectrum of 
the generated modes for different magnetic field spectra. 

4-1 Time Evolution 

Taking the divergence and time derivative of eq. (18), and using eqs. (12), (13), (19) 
and (20), we obtain an equation for the evolution of the divergence of the velocity, 

,.{«,,«V.v))-lf^2v.v=^. (51) 

where 

Q = [V X V X (v X B//o)] X B//0 + (V x B//o) x [V x (v x B//o)] , (52) 

and Bjyo is today's value for the force-free component of the magnetic field. Since the 
perturbation scheme assumes |B//o| » |b|, |B//o| ^ |Bo|. To leading order in this 
perturbation scheme, the time evolution is determined by jBy^joP with no contribution 
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from b, while the spectral dependence is linear in b. When deriving the time evolution 
and the magnetic Jeans length in this section, we may use the fact that |B/y^o| ~ |Bo| to 
omit the // label, but when deriving the spectrum in §4.2, it is important to keep the 
distinction. 

Substituting the explicit time dependence given by eqs. (6) and (7), we can re-write 
eq. (51) in the form 

v.(a„v.|.v-i,v).(^)^Z_S. 

The evolution equation for each velocity mode v(k) is then obtained by inserting Fourier 
expressions for the velocity field (as in cq. (27)), and for the force-free background field in 
analogy to eq. (28). Calling Bj-j-(k) the Fourier transform of Bj-j-o(x, t) as in eq. (28), we 
get: 

k-(t''du + 2tdt-^\ HKt)^ (54) 



P J d^kid^ka F[k,ki,k2,%(ki),%(k2),v(k-ki-k2,t)], 



where 



^2 

f3^ '° 



47rpo 247r2poG' ' 
F = -{k\B, ■ B2){k - k2)i - k\\i - k2) • Bi]S2i 

+ 2(k • B2)[(k • Bi)A;2i - (k2 • Bi)/c,]} {;,(k - ki - k2, t) , 

Bi = B//(ki), and B2 = B//(k2) . For notational convenience, we define a logarithmic 
time variable T = Int, and further define the operator Gi such that 



j Sk^Sk2F = Gi{k:\)v,{\,T), 



I.e., 

Gi = 



= - j (^'k^d^k2dh{e(B^ ■ B2)(fc - k2)^ - fc'[(k - k2) ■ Bi]S2i 

+ 2(k • B2)[(k • Bi)A;2i - (k2 • B^)h]} 5{k - ki - k2 - 1) . 
Upon using this definition, eq. (19) then reads 



dTT + dr-^ i5||(k,T) = ^G,(k:l)i},(l,r). (55) 
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Since the ensemble average of both sides of eq. (55) vanish, we compute the evolution of 
the quadratic quantities. At the expense of yet further algebra, one obtains the following 
three equations: 



dT{\v\l (k, T)|2) = 2Re{v\i (k, T)dTvlik, T)) , (56a) 
dTRe{v^^{k,T)dTmk,T)) = l(|^||(k,T)p) -i?e(^||(k,T)aT^i(k,T)) 

+(|aT5||(k,T)p) + ^Re{G,{k : 1)5,(1, T){5^(k, T)) , (566) 

dT{\dTV\\{KT)\') = ^i?e(i;||(k,r)aT5y(k,T)) -2(|aT5||(k,T)|2) 

+ ^i?e(G,(k : l)v,{l,T)dTvl{k,T)) , (56c) 

where Re stands for the real part. These equations look more complex than they really 
are because the recurring term (Gvv*) can be calculated easily if one notes that {Bv{T)) 
and J G{v\\v±) vanish (see Appendix B). Thus, 

(^.(k : \)v,{\)vl{k)) = -l{B}fo)k'{\mk)\') , 

and 

(G,(k : l)d,{\)dTmk)) = ~{Bjf,)e{v^\{k)dTq{k)) , 

where 

(Bjfo) = 47r J dk k^B)f{k) . 
These considerations allow us to write a single evolution equation for (|i;||(k, T)p), 

dTTT + ^Ott + + \ae\dT + (\ak^ - ^"j |(|5||(k,r)|2) =0, (57) 



where 

' ^'ffo 



a 



(^ffo) 



In addition, we have the initial conditions (at T = Tree): 

(i)y(k,T = Trec)) =0, 

5T(|^;||(k,r)|2)|^_ = 2i?e({;||(k,Trec)aTi;i(k,Tree)) = 0, 

dTT{\v\\{k,T)\^)\^^^^^^ = 2(|aTi;||(k,Tree)|') = 2[trec^^||rms(/i;,irec)]', 
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where i) = dfV and v\\T.^^{k^ t) = •y/(|i;||(k, These follow from the assumptfon that the 
velocity is zero everywhere at t = tree- In the last initial condition above, 'C||rms(^) ^rec) will 
be obtained later upon using eq. (18). 

Returning to our original time variable t = exp(T), and again defining time derivatives 
with respect to the ordinary time t, the solution for the rms velocity can be written as 
follows: 



(i) For k^kB : 



'^^||rms(^; — '^||rms(^) ^rec) ( 7 

\ ''rec 



(ii) For k = ks 



(58a) 



,{k,t) = vn 



t \~2 / t 

i(^) ^rec)^rec ( ~ ] In I ~ 



(586) 



where m(fc) = \/25 — 12aP = 5^/l — (k/kBY- The transition from stable {w? < 0) to 
unstable (m^ > 0) modes occurs at A; = A;^, and we therefore find the magnetic Jeans 
length to be 



27r „ /12a _ 2 (5?/o)'/' 



ks V 25 



5 .„VG ■ '''' 
Modes with length scales greater than \b {k < kB) are unstable, while those on scales 
smaller than {k > kB) undergo damped oscillations. As in Peebles (1980, 1993), one 
can roughly estimate the magnetic Jeans length by using the Alfven speed, va = B/yJAirp, 
in place of the sound speed in the expression for the ordinary Jeans length, which gives 



B 1 



very close to the expression derived in eq. (59) . 
We can now solve for 



5rms(A:,i) = V(l^(k,OP) 
by using the corresponding solution ■C||rms(^) ^)) with the initial condition 



'^rms(^) ^rec) — ) 
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and eq. (12). The detailed derivation of 5rms(^)^) cis a function of Vrmsikjt) is provided in 
Appendix C with the result 

^^Ormsl/e,tj-/e ^^^^ . 



Thus, we can write the solution, 



where, for k ks, 



T{k,t)= 2 
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\ ^rec / 


m 





(m+l) 
6 



for k — ksi 



and, again. 



T(/c,t) = 36( 

Irec 



C{k) 



1 - - In 
6 



rec ~ 
V\\ 



i ^rec) ^ 



(60a) 



(606) 



(61) 



In the form for 5rms(^5 1) defined above, r(A;, t) can be thought of as a transfer function 
which evolves an initial spectrum c{k) from recombination to a later time t. We now discuss 
the behavior of T(fc, t) and leave a discussion of c(/c), which depends on the particular form 
of the magnetic field spectrum, to the next section (§4.2). 

An interesting feature of the solution is that T{k, t) is independent of the spectrum 
of the magnetic field (as long as its statistics are Gaussian) and only depends on the 
combinations k/ks and t/trec- Therefore, the solution can be easily rescaled for different 
choices of magnetic field strength. The time evolution reduces to the usual t^^^ growing 
and decaying modes in the limit k/ks 0. In this limit, m ~ 5 and the growth of 
perturbations is nearly independent of k, so that the final spectrum is mostly proportional 
to the initial spectrum given by c{k) as in §3. As k ^ ks from below, the growth of 
perturbations decreases and, at k = ks, the solution changes from unstable (for k < ks) 
to damped oscillatory (for k > ks)- 

In principle, the time evolution derived above is only valid for modes which correspond 
to scales smaller than the Hubble radius at recombination, k > k^ec — 27r/100 Mpc. Modes 
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on comoving scales between the Hubble radius at recombination and the Hubble radius 
today have similar but delayed time evolution, since these modes start growing after they 
enter the Hubble radius. As we show below, the spectrum of perturbations, c{k), is fairly 
steep for small k, so that the power for small k modes is negligible, and we can approximate 
the power for k < k^ec to be zero. Alternatively, we can estimate the effect of the delayed 
evolution by using t/tenter(fc) instead of t/trec in ^(k^t) for modes with k < kj-ec, where 
tenterik) is the time a /c-mode enters the Hubble radius, tenter(^) = trecikrec/k)^. 

In Fig. 1, we plot T{k,t/trec) for k > kj-ec and t(/c, t/tenter) for k < krec at different 
redshifts. (The sharp discontinuity at A;rec is an artifact of the approximation that recom- 
bination happened instantaneously.) We can see that, independent of the spectrum of the 
magnetic field, the power on scales k < k^ec is suppressed (by the delayed growth) as is 
the power on scales k > ks- Therefore, while magnetic fields do not generate significant 
clustering on scales larger than ~ 27rfc~J and smaller than ~ 27r/i;^^, they can have signif- 
icant influence on the formation of structure between these two scales, depending on the 
strength and spectrum of magnetic fields at recombination. 

If the universe were baryon-dominated with infinite conductivity, magnetic fields could 
deter the growth of perturbations on scales smaller than As; hence, a constraint on the 
strength of magnetic fields could be derived by the observation that structures do form on 
a given scale above the ordinary Jeans length Aj. For example, for galaxies to form in a 
baryon-dominated, infinitely conducting fiat universe, we would require Xb ^ Ig — 1 Mpc 
and, therefore, (Bq)^/"^ ^ 4 x 10~® G. However, this result does not take into account the 
presence of neutral hydrogen or non-baryonic dark matter. When these components are 
present, modes on scales below Xb are also unstable. 

Finally, the time evolution of incompressible modes can be derived in an analogous 
manner to the derivation of the compressible mode evolution. We should also note that the 
incompressible modes do not affect the compressible modes just studied, and vice-versa, 
an assertion which can be demonstrated by considering eq. (20) in detail. 



In order to study the incompressible modes, we take the curl and the time derivative 
of eq. (18), use eqs. (19)-(20), and find 



which is the incompressible analog of eq. (53). We then take Fourier transforms as in eqs. 




(62) 
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(27)-(28), and obtain the analog of eq. (57) for the incompressible mode: 



dTTT + 2,dTT + (f + \c^k^^dT + i^-ae + ^) }(l^±(k, m = , (63) 

where kv^Ck^trec) = |k x v(k,ti.ec)|- 

Since the calculation is essentially the same as that described above for the compress- 
ible case, we only present the solution for the incompressible case here: 

(i) for ak'^ 1/6: 



{k,t) 

— ^_Lrms (^) ^rec) I j. 



rec 



t V f J 

tr 



(64a) 



(ii) for ak'^ — 1/6: 



^_Lrms(^, t) = V±rms{k, trec)^rec ( 7^ ) J—j : (64&) 

V'-rec/ V'-rec/ 

where v±rms{k,t) = y/{\v±{k, t)\^), p{k) = Vl - 6aA;2/6, and a = {BjfQ)/127r'^plG is 
the same as in the compressible case. Unlike the compressible modes, the incompressible 
modes have no growing solution, only decaying or damped oscillatory solutions. 

4-2 The Spectrum of Compressible Modes 

In order to obtain the spectral dependence of the generated compressible modes, we 
need to find 5||rms(^) ^rec) or c{k) for a given magnetic field spectrum at recombination. 
From eq. (18): 

^ . . , , V-{[(VX B// rec) X brec] + [(V X brec) X B//rec]} 



47rp5(trec)-Rrec 

V • {[(V X Bffo) X bp] + [(V X bo) X B^^o]} 

47rpo^rec 



(65) 



and we need to specify spectral dependencies of both B/jo and bg. 

In what follows, we discuss two Ansatze for the functional form of the magnetic field 
spectrum. We study the solutions for the cases in which the perturbations about the force- 
free field are either a delta function or power laws. In both cases, we choose a delta function 
for the background force-free spectrum to make the problem tractable analytically. This 
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choice adequately describes the force-free background as long as the scale of perturbations 
onto this background are much smaller than the background scale. 

Assuming the force-free background magnetic field spectrum has a spectrum peaked 
at some scale, Aq = Stt/Zcq, we can describe it by 

Bh(k) ^ (Bho)^-^^ ■ (66) 

In Appendix D, we show that for this choice of B'j^{k) the initial acceleration can be 
written as: 

+ + 2k^k^ + 16k^kl + 2k^)kl + kl- 2k'^kl + 2k^kl - k^\ , (67) 



and ^{ki) is specified in the next two sections. 

4..2.1 Delta Function Spectra for h^{k) 

If we assume that the deviations from the force- free background are also peaked around 
a scale kh, then we can write 

(68) 

where 

{hi) =At^ j dkk%^{k) . 
The spectrum of generated compressible modes is then given by 



^Jn^hM) 2kh 



^5||rms(^,W) = — 52 7^=^VFs{k,ko,kb) (69) 

(47r)'^P0-Rrec Vkbko 

for kb — ko < k < kb + ko, where Fs{k, ko, kt) is given in Appendix D, eq. (D7). Using eqs. 
(59) and (61), we find for the "initial" density spectrum 



Although kf) ^ ko within our assumption of scale separation, we may take the limit 
kb ko to make a connection with §3.2.1. In this limit, we recover eq. (37) with an overall 
factor of 2 difference due to the definitions of (B^^q), (6q), and (Bq). 
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Given the linear relation between the generated compressible modes and the magnetic 
field perturbations about the force-free background, (i.e., 5 oc 6), it is not surprising that 
the range of excited modes of 5 is very narrow {kb — ko < k < ki, + ko) for the case of a 
delta function spectrum for b. In this scenario only structure with k ^kb may form. 

Although this scenario is not very realistic, we briefly discuss the formation of non- 
linear structure for the following limiting cases: kb — fcs and kb ^ ks- For kb — ks, we 
set A{kB,tni) — c{kB)T{kB, tni)^J^2T:'^k^ /V = 1 and find that tni satisfies 

- (l-,,) = l + -ln - , (71) 



where 75 = 3.8 x ^k^/kB^ ^-^ffo) /(^o)- Since tni tree as 75 — >• 0, the smaller 
ko/kB the earlier non- linear structure on scale kB can form. 

As kb/kB decreases from unity, the cutoff density spectrum becomes kb + ko instead 
of kB- For kb <^ kB, can approximate m ~ 5 in eq. (60a) and write: 



2 • ■ ^ ' ■2//r2\\5 



.hj \kbj y {%) 

The time dependence in this case is similar to the examples discussed in §3.2.1. 

The narrow range of excited modes limits the applicability of these results to physically 
interesting scenarios. We therefore move on to the more physically relevant case of power 
law spectra for b{k). 

4-2.2 Power Law Spectra for h^{k) 
The more general case of a power law spectrum of magnetic field perturbations can 
be studied if we assume that, for Ajmin < ^ < ^max, 

h^{k) = Ck''. (73) 

where C = (9 + 3) {hi) (A;^+3 - C+n)"'- For k^r. « /^max, C ~ (5 + 3) {hi) k'f'^ In 
Appendix D, we show that for this Ansatz (eq. (D8)) 

k'iUsik, W) ^ {Bjfo) {hi) ^ + 3) F,{k, ko, /c^ax) (74) 

where Fq{k,ko,kma.x.) is given by eq. (D9). Therefore, we find that 
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3 



and 

96 V 2 \ ko J kg 



where 



In Fig. 2, we plot Fq{k, ko, /cmax) for different choices of q with fixed /co/^max = 10~^. 
If we take the limit A;o/^max — > in eq. (D5) (note that the limits of integration depend 
on ko), Fq can be approximated by 



Fg ~ 16 



ko V f k ^^+' 



k \ k 

"■max / \ "^max 



as it is clearly shown in Fig. 2. We can, therefore, write 



and 



9 + 7 
2 



t) - S ( ) 1^1 r(k, t) . (78) 



J \ ""max , 

The spectrum of generated density perturbations in this case depends both on the 
magnetic field spectral index {P^ oc k^^"^, for /c <^ /cg) and the ultraviolet cutoff /Cmax- The 
amplitude of density perturbations is determined by ks, fcmax, and r, the ratio between 
the energy density in magnetic perturbations relative to that of the background field. We 
have concentrated on > — 3 such that the dependence on /cmin is negligible. The results 
can be easily extended for cases in which A;inin plays an important role. 

As an example, we choose the case of = and plot Ao(A;, t)/^Jr at different redshifts 
in Fig. 3, asssuming ks = Aj^ax, ko = IQ-^kme,^, and A;„iin = 10~^A;max- For k < k^e,^, 
Pq OC k^ and the peak power is at A; ~ Aj^ax = ks- In this example, the behavior at small k 
is too steep to fit large scale structure observations and structure will be formed primarily 
around the scale A;max or kB- To broaden the infiuence of magnetic fields on the range of 
linear structure a spectral index q ~ —3 is necessary (see Ratra (1992a,b), for an example 
of g = — 3 primordial magnetic field spectrum). 

We now discuss some examples within the large parameter space (A;inax) ko, ks, and 
r) that are relevant to the formation of structure in the universe. As long as q > —3, the 
first scales to become non-linear have k ~ min(A;inax) A;^). We again define tni{k), such 
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that Aq{k,tni) = 1. For A;inax + A;o > ks, ks is the largest mode to be excited. Setting 
^qi^Bitnl) = 1 also leads to eq. (71) but with 



3.8 X 10-2 //C^ax\'^' 



7<z 



^/r(gT3) V 



Since as 7g 0, t-ni —>■ tree, for /Cmax + ko > ks the closer A;inax is to kB, the earlier 
structures on scales ks can form. 

To simplify our discussion, we set g = and r = 1 and vary /cmax, ^s, and ko (different 
choices of q and r can be accommodated by making the appropriate rescaling of /cmax)- 
We find that for /cmax — ^s, the first non- linear objects could form as early as redshift 
Znl{kB) — 300 and the scale of non-linearity today would be kni{to) ~ 0.15 kB- If we set 
knl{to) = 27r/8 Mpc, then fcmax ^ kB — 27r/l Mpc, which satisfies eq. (10). Since the 
background field is a delta function the appropriate bound on the parameters due to the 
upper limit on magnetic fields on Mpc scales is given by eq. (43) with ko substituted for 
k^ (as long as r <S 1). In this example, eq. (43) is satisfied if /cq = 27r/2 Mpc. If instead, 
we choose kni{to) — 2tt/2 Mpc, then /Cmax — ks — 27r/0.3 Mpc which satisfies eqs. (10) 
and (43) with {B^y/^ ^ 10"^ G and B{lg,to) ~ 10'^^ G (for ko ~ 27r/2 Mpc). 

The time at which modes with k = ks become non-linear, tniiks), increases as 
kmax/ks grows from unity. For magnetic fields to play a role in structure formation, 
tnii^B) must be less than the age of the universe, to, which implies that /cmax ~ 6A;b. 

As kma.x/kB decreases, so does tni{kB) up to Ajmax + ^o — ^b, where the cutoff changes 
from kB to A;niax + ^o- For fcmax + < ^S) we follow tn/(^max) using our solution (eq. 
(60a)) in the limit k <^ kB (m ~ 5) which implies: 

tni{k)Y 1.5 / kB Y k 

max 



tree J ^Jr{q-\-^) ^ \ "'max / \ 

From eq. (79), requiring tnz(^max) ~ leads to the constraint fcmax ~ lQ~^kB- This 
estimate helps define the range for which magnetic fields can make non-linear structures, 
i.e., 10-2/j^ < k^.^^ ^ Qks- For example, we can choose /cmax — O.l/cs, which gives 
^nz(/Jmax) — 7. Choosiug /cmax — 27r/0.8 Mpc, then ks — 27r/80 kpc and all the observa- 
tional constraints are satisfied, with (S^)^/^ ^ 3 x lO-^ G and B{lg,to) ~ 10"^° G. On 
the other hand, if kma.x — lO~'^kB — 27r/8 Mpc, then tni{kma.x) — to and kB = 27r/0.08 
Mpc which satisfies eq. (10) but violates eq. (49). 
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We see that with a more reahstic choice of power spectra for magnetic fields at recom- 
bination, non-hnear structures in a cosmologically relevant range of scales can be formed 
at reasonable redshifts for field strengths below the observational upper limits. Although 
the number of parameters in this discussion is large, the physically relevant range is quite 
limited. As we improve our understanding of the evolution of magnetic fields from the early 
universe up to the epoch of galaxy formation, the parameter space will become more con- 
strained. For instance, in the absence of primordial vorticity, magnetic field perturbations 
are damped up to the Silk damping scale and /cmax ~ 27r/10 Mpc (Jedamzik, Katalinic, & 
Olinto 1996). 

5. Conclusions 

We have studied the effects of magnetic fields present at recombination on the origin 
and evolution of density perturbations and peculiar velocities. We find there are generic 
features of the generated density perturbations which are largely independent of the as- 
sumed spectrum of the primordial magnetic field. The first conclusion we can draw is 
that, within our linear perturbation approach, magnetic fields cannot explain the observed 
galaxy power spectrum on large scales, since the generated spectrum is too steep for small 
k. To reproduce the observed P oc k spectrum on large scales, we need b^{k) oc for 
small k with amplitudes that violate the observational constraints on (Bq). 

Another generic feature is the cutoff introduced by the magnetic Jeans length. This 
cutoff limits the amplitude of the power spectrum for any choice of magnetic field strength. 
As the magnetic field strength increases, the amplitude for a given density perturbation 
mode rises, but, simultaneously, the magnetic Jeans cutoff moves to smaller k. The net 
effect is that the peak amplitude for the resulting density power spectrum does not neces- 
sarily increase as one increases (Bq). 

Since the generated spectrum falls sharply at small k and is cut off at large /c, magnetic 
fields generally produce a peak in the density spectrum over a narrow range of wavenum- 
bers. For this peak to be of relevance to the formation of structure, the amplitude A{k, to) 
must be ^ 1 for scales 27r/c~^ ^ 8 Mpc. The peak amplitude is sensitive to the assumed 
spectrum for the primordial magnetic field; the smaller /cmax, the stronger the variance 
in the density perturbations. In particular, if fcmax ^ ks, density variances well above 
unity can be obtained with relatively small magnetic fields. Depending on A;^ and A;inax5 
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objects from galaxy scales down to a first generation of massive stars can be formed. As 
the variance reaches unity our calculations break down, and another conclusion can be 
drawn: when magnetic fields are important, their effects are mostly non- linear in nature. 
Our linear calculations can, however, be used to estimate the epoch of non-linear collapse 
of different mass scales. 

In the present work, we have focused on a purely baryonic universe. In a subsequent 
paper, we study the evolution of density perturbations when non-baryonic dark matter is 
the dominant component of the universe. When we include non-relativistic (cold) dark 
matter, we find that the amplitude of the density perturbations decreases as fibaryon de- 
creases for a fixed magnetic field strength. On the other hand, the perturbations become 
unstable for all wavenumbers: since cold dark matter does not couple to the magnetic 
fields, no damped oscillatory modes can be sustained. In this case, the perturbation am- 
plitude A{k) flattens out for large k rather than being cutoff at fcg. Baryons will still show 
some resistance to clumping on small scales due to the magnetic fleld; this may segregate 
baryons from dark matter, introducing a source of bias on small scales. 

In the case of hot dark matter the growth of baryon perturbations on small scales is 
slower due to neutrino free-streaming. Therefore, the initial perturbations need to be much 
larger for the final variance to be greater than 1. Clearly, the hot dark matter scenario 
would greatly benefit from a peak in the variance at large k so that structures can form 
on scales smaller than the neutrino free-streaming length and populate the problematic 
empty voids. 

Finally, we conclude by noting that magnetic fields most likely play a dynamical role 
in the formation of galaxies and clusters of galaxies even if the original perturbations were 
caused by other sources. 
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Figure Captions 

Fig. 1: The transfer function T{k,t) plotted versus k/kB for redshifts z = 0,1,5,10,100, 
and 1000. We plot r(/c,t/trec) for k > krec and t(A;, t/tenter) for k < /crec with tenter{k) = 

tvec{krec/k)^. (We Set krec = 0.01/cb.) 

Fig. 2: The function Fq for the power law Ansatz as a function of k/kmax, for ko = 

10 ^/Smax) ^min = 10 ^^max) and q = —3, —2, —1, 0, 1, 2. 

Fig. 3: The variance for the power law Ansatz with q = 0, Aq/a/t, as a function of k/ks 
for redshifts -2 = 0, 1, 5, 10, 100, and 1000. (We set /cmax = lO^/cmin = lO^/co as in Fig. 2.) 
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Appendix A 

In this Appendix, we derive eq. (34), apply them to the two Ansatze for B(k), and 
derive eqs. (36) and (45). We begin with eq. (31), 

V • V(X, tree) = «V • [(V X Bq) X Bq] , (^1) 

where a = l/47rpo-Rrec- By using the Fourier expressions for v and Bq, we obtain 

ik ■ ^(k, W) = « y d^^i k • {ki [B(ki) ■ B(k - ki)] - B(ki)[ki • B(k - ki)]} , {A2) 

where B = B throughout this Appendix. If we choose our coordinate system in such a 
way that k is along the 2;-axis, 



d^ki =4:77 dkikj / dfx , 



where = cos 9 with the angle between ki and the z-axis. The integration range in eq. 
(A2) should be taken carefully since the integrand had 5{k — ki — k2) before integrated 
over k2 with the condition of A;min < < ^max for z = 1,2. In other words, there is a 
constraint on the angle to be integrated over, depending on the magnitude of ki. By taking 
this constraint into account, the integration range of eq. (A2) depends on the magnitude 
of k, and can be shown to be 

(1) For < /c < kmin. 



dkid/j, — / dki / d/j, + dki / d/j, 

/•femax /•! 

+ A /^V^!±^i^^^- ^^^^ 



2fcfci 



(2) For /cmin <k<k^ 
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/pk-kmin rl rk + kjnin r 2fefei 

dkid/j, — / dki / d/j, + dki / d/j, 

^ Lk ^^V/^ + A V^!±^iz^ • ^^^^ 

-'fc + femin -^-l -'Kmax-K jtlj^^ 

(3) For A;niax ^ k ^ ^max + ^min" 

fvinin nrnn •'K — Kmin oTTin 

(4) For /Cmax + ^min ^ ^ ^ 2/Cmax- 

dkidn = / d/ci /,2^,2_,2 ci/i • (^6) 

^min 



%1 Lift/ci 



If we take the ensemble average of both sides of eq. (A2), and use eq. (29), we obtain 

(k-i^(k,trec)) =0. 

To calculate 

^rms(^; ^rec); take the square of both sides of eq. (A2) and obtain 
|k ■ ^{k, = y d^kid^k2{{-k ■ B*(ki))[ki ■ B*(k - ki)] + (k ■ ki)[B*(ki) ■ B*(k - ki)]} 

X {(-k • B(k2))[k2 • B(k - k2)] + (k • k2)[B(k2) • B(k - k2)]}, {A7) 



We then take the ensemble average of both sides of eq. (A7), and express the fourth-order 
correlation of S as a sum of a product of second order correlations by assuming a Gaussian 
statistics for B. Since each fourth-order correlation gives three products of second-order 
ones, we obtain a total of twelve terms on the right hand side of eq. (A7). Applying eq. 
(29) to these twelve terms yields terms proportional to 5(k), 5(ki — k2), and 5(k — ki — k2), 
respectively. Four of these terms (which contain 5{k)) vanish upon integration and they 
will not be written here; retaining the remaining eight terms leads to the expression 
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(|k • ir(k, t,ec)|') = y d^k^d^k^iiik ■ B*(ki)k • B(k2))(ki • B*(k - ki)k2 • B(k - k2)) 

+ (k • ki)(k . k2)(s;(ki)s,(k2))(s*(k - ki)s,(k - k2)) 

- (k • k2)(k • B*(ki)Si(k2))(ki • B*(k - ki)Si(k - ks)) 

- (k • ki)(k • B(k2)S*(ki))(k2 • B(k - k2)5*(k - ki)) 

+ ((k-B*(ki)k2 ■B(k-k2))(k-B(k2)ki ■B*(k-ki)) 

+ (k • ki)(k • k2)(s;(ki)s,(k - k2))(sr(k - ki)s,(k2)) 

- (k • k2)(k • B*(ki)S,(k - k2))(ki • B*(k - ki)S,(k2)) 

- (k • ki)(k • B(k2)S*(k - ki))(k2 • B(k - k2)S*(ki))} 

= ^^'f^||rms(k,trec) , (^8) 

where k • v(k, tree) = (k, tree)- The integral of these eight terms is given by (using eq. 
(29)), 

,2^2 „ , a^VP „ ,2 , S2(A;i)S2(|k-ki|) 
^ ^||rms(^>W) = / dkik^ dfl 

\ / J krain */ 1 I 1 I 

X {2k^kiiJ, + k'^klil - 5//^) + 2kkln^] , {A9) 

where the identity 6{k = 0) = V/{27t)^ has been used. Eqs. (A8) and (A9) give eq. (34). 
For the delta function Ansatz, B'^{k) — {BQ)6{k — k^)/AT:k'^, and therefore 

Eq. (36) follows from this result. 

• 2 ~ 

We now evaluate Vi^^^.^^ for the case of a power law spectrum for B'^{k), as in eq. (44). 
For simplicity, we assume /cmin = 0. The resulting integral can be solved analytically if we 
use /c//cmax a small parameter. In such an expansion, the leading order term turns out 
to be (k/kiaax)^ for Q = —2, and (A;/A;niax)^ for q = —1, 0, 1, 2, 3, 4, and 6. Finally, we 
find the following solution for the q > —1 case: 
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"'max 



^^■^||rms(^' ^rec) 



{Any 



15(2g + 3) 



l + O 



(All) 



Therefore, 



|k- V(k,trec)|'') = k'^vZmsik.trec) 



/o2^2 22{q + 3) 



{A12) 



where we used (B^) = Ati j'^^Z ^^^^ - ^^mtx/(? + 3); this leads directly to eq. (45). 

Appendix B 

In this Appendix, we demonstrate that (i?i.ec'f^|| (^)) = {BrecV±{t)) — 0. We begin with 
the evolution equation for v^\, obtained through eq. (54) 

k(^t^dtt + 2tdt-^^ {)||(k,t)=/3 J (i3kid3k2F(k,ki,k2,B(ki),B(k2),v(k-ki-k2,t)) 

m 

With the use of the Green function for this equation, 

3e{t - 1') 



G{t, t') = 



5kt' 



_1 4 



{B2) 



(where 9 is the step function), one can solve eq. (Bl) formally, given the initial conditions 

V\\{k,trec) = 

and 

dtV\\{k,t)\t^trec = ^\\{k,trec) ■ 

This Green's function solution is given by 

Strec'^ll {k, ^rec) 



V[\{k,t) = 



— 1 4 



-rec 

T 



^rec 

T 



+ 



(3 dt' d^kid^k2G{t, t')F{k, ki, k2, B(ki), B(k2), v(k - ki - k2, t')) (S3) 

tree 

This equation can be solved by iteration, yielding a Born-series solution, for which the 
zeroth-order iterate v?.{k,t) is given by 



3trec'^|| {k, tree) 



rec 

T 



1 4 
3 // \ 3 



(S4) 
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We first note that {v'^^{k,t)Bj.ec) = 0. This result is obtained simply by noting that 
V\\{k,trec) is quadratic in Srec (cf. eq. (55)), so that v^^{k,t)BT.f;(. is an odd product of Sj-ec. 
Hence, as long as By-^c is a random variable governed by Gaussian statistics, ({Sy (/c, t)i?rec) 
must vanish since all odd moments of a Gaussian random variable vanish. 

It is now straightforward to demonstrate that {v\\ (k, t)B,-ec) = to all orders by simply 
computing the next-order iterates: At the n*^ iteration stage, v^{k,t) is the sum of even 
products of Srec (because F is a quadratic form in B^ec)', hence, v^{k,t)Bj-ec must be 
always an odd product of -Bj-ec, and therefore {v'!^{k,t)Brec) = 0, as before. Thus, under 
the assumption that the Born series expansion converges, we obtain a portion of our desired 
result, namely that 

({;||(A;,t)5,ec) = 0. 

It is now readily shown that the same result obtains for the incompressible flow, v±, 
by simply repeating the above calculation, but now projecting out the incompressible 
component. Thus, we also obtain {v±{k,t)Bj.f;c) = 0, and hence we have the ultimately 
desired result that 

{v{k,t)B,^,) = 0. 
Appendix C 

In this Appendix, we show how to obtain drnis{k,t) from 'Urms(^)^) in the force-free 
case of §4. We start with eq. (55) in the main text, which reads 

(a^^ + dT- ^)5||(k,T) = ^G,(k : 1)5,(1, T) . (CI) 
Upon using T = In t, this is equivalent to 

{t%t + 2tdt - ^)5||(k,t) = ^Gi{k : 1)5,(1, t) . (C2) 

The rms solution to (C2) was shown to be given by eq. (58) in the main text, which we 
rewrite as 

^'||,rms(^,^) = 5||,rms(^,^rec)l^(^) , (C3) 

where 
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Vrms{k, t) = •y/<~[v|i'(k7^JP3^, and m = -\/25 — 12ak'^. 

First, we multiply both sides of eq. (C2) by v*^{'k,t'), with t' ^ t, and then take an 
average to get the following equation: 

{t'du + 2tdt - ^)(^||(k,t)i)|(k,t')) = ^(Gi(k : l)v,{\,t)vl{k,t')) . {C4) 

Now, after repeating the same calculation that was used to get eq. (57) from eq. (56), it 
can be shown that 

(G,(k : l)vi{l,t)vl{h,t')) = -|(S2)A;3(^||(k,t){i|(k,t')) • (C5) 
Using eq. (C5), we readily obtain the solution to eq. (C4), namely 

{mk,t)vl{k,t')) = {^{k,t,,,)v;{k,t'))V{t) . (C6) 

Next, let us take the complex conjugate of both sides of eq. (C2), multiply by v\\(k.,trec)i 
and then take an average; that gives 

(5ff(k,t)5||(k,W)) = (|5||(k,W)|')^(t) . (C7) 

From eqs. (C3) and (C7), 

(^y(k,t)5||(k,trec)) = y^(|£i||(k,W)P)V^(|5||(k,t)P), 

or 

(^^(k,t0^ii(k, w)) = 7(|^||(k,w)P)7(hi(k,t')P)- (cs) 

Upon using eq. (C8) in eq. (C6) 

(^ll(k,t)^|(k,t')) = ^(|^i||(k,W)P)7(h(k,t')P)V^W- {C9) 
Finally, with the help of eq. (C3), eq. (C9) becomes 



(i;||(k,t)i;|f(k,0) = ^(h(k,t)P)^(h(k,tOI^). (^^10) 
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On the other hand, eq. (12) in the main text gives us the following equation: 

By integrating eq. (Cll), the solution for the density fluctuations is 

~6{k,t)^-ik f (C12) 



^re 



Rit') 



By multiplying by the complex conjugate of eq. (C12), and then taking its average, we 
obtain 



/•t Jj. rt ^j. 



(C13) 



where eq. (CIO) has been used; taking the square root of eq. (C13), we obtain 

7(|5(k,t)P) = /c^' ^^(|^ll(k,ti)P). (CM) 



In other words, 

/* dt 
-^^Vrm8{k,ti) . (C15) 

The differential equation for ^rms follows via taking the time derivative of eq. (C15), 

Appendix D 

In this Appendix, we show how to derive eqs. (69) and (74). We begin with eq. (65), 
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V • v(x,trec) = a{V • [(V X B//o) X bo] + V • [(V X bo) X B//o]} , (Dl) 



where, again, a = l/47rpo-R?ec- By using the Fourier expressions for v, B//05 and bo we 
obtain 

ik ■ ^k, tree) = aj d^k^ {A;'[b(ki) • B//(k - ki)] - [(k + ki) • B//(k - ki)][k • b(ki)] } , 

m 

where B//(k) is the Fourier transform of B//o and b(k) is the transform of bo- The 

ensemble averages of both sides of eq. (D2) vanish. Again, we compute 'i'rms(^, ^rec) by 
taking the square of both sides of eq. (D2) and then ensemble average, with the following 
result: 

'k ■ J-1 k - ki M 



|k-v(k,trec)|^) / dkikl diJi 

*^ fcmin *^ 1 



X [k^ + k'^klil -11'^ + - k^kiii{l + At^)} . {D3) 

The polynomial integrand in eq. (D3) is not the same as in eq. (A9) since B// and b are 
assumed to be uncorrelated and separated in A;-space (scale separation). This assumption 
changes the non-trivial contributions in eq. (A8). 

For the force-free background magnetic field, we assume a delta function spectrum 
and write 

BhW^iBj,,}^-^^. (DA) 

where we assumed that ko is small compared to the wavenumber of the fluctuations about 
the force-free background. Eq. (D3) now takes the form 

(|k • Hk. W)P) = ^^0^ fj^^ dk,^^ i^kl - {Akl + 2k')kt + {6k', + 2k\l)kt 
+ + 2k'^kl + IQk'^kl + 2k^)kl + kl- 2k^k^ + 2k^kl - /c^j . (D5) 

For a delta function Ansatz for lP{k), 
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only modes satisfying kb — ko < k < kb + ko are excited. With the use of eq. (D5), the 
spectrum of the generated compressible modes is given by 

• 2 4 a^Vk^ 

fe'^llrms(fe.W) ^ (4^)4fc^^fc^ {Bjfo){bl)Fs{k,ko,kb), {D6) 

where 

w (h I. ^ — 1 I (^b ~ ^o)'^ ~ ^^^(^b ~ ^o)(^6 ~ ^o) + ^^^(^b + ^o) ~ (r>7\ 

from which eq. (69) follows. 

In the limit kb ko and for k <S 2kb, we recover the result in eq. (AlO) with an 
overall factor of 4 difference due to the relation between (S^^q), (69), and (Bq). 

The power-law Ansatz can be written as 



b^{ki) = Ckl 



q 

1 



for k^^ <ki< /c^ax, and (bl) = AttC - + 3) ~ AttC k^+,l/{q + 3) for 

^min ^ ^max- Substituting this Ansatz in eq. (D5), the integration can be performed 
within the following range: for k > k^g.^ — ko, the upper limit on ki is max(A;i) = /^max; 
and for k < /cmax — ko, max(A;i) = k + ko; while for k < \/^o~+~^Iiin' lower limit is 
min(A;i) = /Cmin, and finally, for k > a/Zcq + k'^^^, min(A;i) = a/A;^ — /cq. The solution for 
specific choices of ko, /cmin, /^max, and q is plotted in Figs. (2) and (3). 



Since ko <^ k^in -C fcmax, the most relevant range of the integration is \/^o~+~^min — 
k < /cmax — ^0- In this range, the integral in eq. (D5) can be written as 

k%msik, tree) ^ i^(^//0)(^0)^ (Q + 3) F,{k, ko, /C^ax) , (DS) 

where 

1 _ 4^1+2^ 64+2^ 
i^g(/e,/eo,/emaxj ^ + 6 + g + 4 

H : J ; {^^) 



q + 2 q 
Ko = /co/^max and K = k/kraax- In the limit kq <^ k, Fq 16 Kq k'"'"^ and, therefore, 

/c^i^llr^^l/c, W) ^ (Sj^o) (^'o) + 3) j . (DIO) 
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This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9412070v2 
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